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ABSTRACT 

Several  researchers  have  recently  developed  new  techniques  that  give  fast 
algorithms  for  the  minimum-cost  flow  problem.  In  this  paper  we  combine 
several  of  these  techniques  to  yield  an  algorithm  running  in 
0(rm  log  log  U  log(nC))  time  on  networks  with  n  vertices,  m  arcs,  maximum 
arc  capacity  U,  and  maximum  arc  cost  magnitude  C.  The  major  techniques  used 
are  the  capacity-scaling  approach  of  Edmonds  and  Karp,  the  excess-scaling 
approach  of  Ahuja  and  Orlin,  the  cost-scaling  approach  of  Goldberg  and  Taijan, 
and  the  dynamic  tree  data  structure  of  Sleator  and  Taijan.  For  nonsparse  graphs 
with  large  maximum  arc  capacity,  we  obtain  a  similar  but  slightly  better  bound. 

We  also  obtain  a  slightly  better  bound  for  the  (uncapacitated)  transportation 
problem.  In  addition,  we  discuss  a  capacity -bounding  approach  to  the 
minimum-cost  flow  problem.  r'' 
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1.  Introduction 

The  minimum-cosi  circulation  problem  calls  for  finding  a  circulation  of  minimum  cost  in  a 
network  whose  arcs  have  flow  capacities  and  costs  per  unit  of  Sow.  Our  framework  for  studying 
this  problem  is  as  follows.  See  e.g.  [9,11,13,18.22].  Let  G  =  (V,E)  be  a  directed  graph  with  ver¬ 
tex  set  V  and  arc  set  E.  We  require  G  to  be  symmetric,  i.e.  (v.w)  e  £  if  and  only  if  (w,v)  e  E.  We 
assume  that  G  contains  no  multiple  arcs;  our  algorithms  easily  extend  to  allow  multiple  arcs. 
Graph  G  is  a  network  if  each  arc  (v,w)  has  a  nonnegative  real-valued  capacity  u(y,w)  and  a  real¬ 
valued  cost  c(v,w).  We  require  that  the  cost  function  be  antisymmetric,  i.e.  c(v,w)  =  -c(*\v)  for 
all  arcs  (v,w)  e  E.  We  denote  by  n,  m,  U,  and  C  the  number  of  vertices,  number  of  arcs,  max¬ 
imum  arc  capacity,  and  maximum  absolute  value  of  an  arc  cost,  respectively.  Time  bounds  con¬ 
taining  U  or  C  are  subject  to  the  assumption  that  all  arc  capacities,  or  all  arc  costs,  respectively, 
are  integral.  For  ease  in  stating  time  bounds,  we  assume  (without  loss  of  generality)  that  Ct2 
and  U  2  4.  All  logarithms  in  this  paper  are  base  two  unless  an  explicit  base  is  given. 

A  pseudoflaw  f  on  a  network  G  is  a  real-valued  function  on  the  arcs  satisfying  the  following 
two  constraints; 

/(v.w)  S  u(v,w)  for  all  (v,w)  e  E  (caj*.  .y  constraint)  (1) 

^  Soan  Mro d  at  Mui|(nnt,  M  L 7..  Cambridge,  MA  1X2139.  Research  partially  Reported  by  an  NSF  Presidential 
Y curt  i  |ivgii|gcr  Feflowrisp,  Co etna  8451517ECS,  grit*  AFORS-88-0088  front  die  Air  Force  Of  See  of  Scientific 
Reward),  end  grants  from  Analog  Devicea.  Apple  Comptaer  Inc.,  and  Prime  Coenpeaer 
Oi.  1— Indian  Institute  of  Technology,  Karpur,  India. 

3  Department  of  Computei  Sdrr.ce,  Stanford  University,  Stanford,  CA  94305.  Reward)  partially  Reported  by  NSF 
Presidential  Young  Investigator  Grant  CCR-8858097  and  ao  IBV  Faculty  Development  Award. 

4  Department  of  Computer  Science,  Princeton  University,  Princeton,  NJ  08544  and  A7ATB-U  Labor»tnrie»,  Murray  Hill, 

NJ  07974  Retea  reft  partially  supported  by  National  Science  Foundation  Grant  DCR-8fi05962  and  Office  of  Naval 
Reward)  Comma  N00014-87-K-0467. 


-2- 


/(v.w)  =  -/ (w,v)  for  all  (v.w)  e  £  (antisymmetry  constraint). 


For  a  pseudoflow  / and  a  vertex  v,  the  balance  bj(y)  at  v,  is  the  net  flow  into  v: 


bj(y)  =  £  f(u,v). 

*K*,*)*  E 


The  cost  of  a  pseudoflow  /is  defined  as  follows: 


cost  if)  —  £  c(v,w)  f(v,w) 

/(».*)>  o 


(2) 


(3) 


(4) 


A  pseudoflow  / is  a  circulation  if  the  following  constraint  is  satisfied: 

bj(y )  =  0  for  every  vertex  v.  (5) 


The  minimum-cost  circulation  problem  is  that  of  finding  a  circulation  of  minimum  cost  in  a 
given  network. 

The  minimum -cost  circulation  problem  has  been  intensively  studied  for  over  thirty  years. 
See  e.g  [9,1 1,1 3,18 ,22].  Among  the  known  algorithms  for  this  problem,  there  are  three  that  have 
the  best  worst-case  time  bounds.  Each  of  these  algorithms  is  best  for  a  different  range  of  the 
parameters  n,  m ,  U,  and  C.  The  algorithms  are  the  0((m  log  UYjm  +  n  logn))-time  method  of 
Edmonds  and  Karp  [6],  the  Olnm  log (nJ/m)  log  (nC))-time  method  of  Goldberg  and  Tatjan  [11], 
and  the  0((m  log«)(m  +/i  logfl))-time  method  of  Oriin  [16].  The  last  of  these  methods  is 
strongly  polynomial  *. 

One  important  idea  is  common  to  all  three  of  these  algorithms,  that  of  scaling  or  successive 
approximation.  Scaling  methods  work  by  solving  a  sequence  of  more-and-more  accurate  approx¬ 
imations  to  the  original  problem.  The  approximations  are  obtained  either  by  relaxing  some  of  the 
numerical  constraints  or  by  ignoring  some  of  the  precision  of  the  numeric  parameters.  Scaling 
was  introduced  by  Edmonds  and  Karp,  whose  algorithm  scales  capacities.  Orlin’s  algorithm  is  a 
refinement  of  that  of  Edmonds  and  Karp  that  combines  capacity  scaling  with  repeated  arc- 

*  A  network  algorithm  it  strongly  polynomial  if  iu  running  time  it  polynomial  in  n  and  m,  assuming  arith¬ 
metic  operation*  take  unit  time,  and  also  polynomial  in  n  ,  m,  k>g(/,  and  krgC  assuming  arithmetic  operations 
lake  lim'  polynomial  in  the  number  of  bits  if  the  operands.  See  121). 
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shrinking.  The  algorithm  of  Goldberg  and  Taijan  scales  costs.  It  relies  crucially  on  the  notion  of 
e-optimality,  introduced  by  Tardos  121]  and  independently  by  Bertsekas  ]3,4]. 

Consideration  of  these  algorithms  suggests  the  question  of  whether  capacity  scaling  and 
cost  scaling  can  be  combined  to  yield  an  algorithm  faster  than  any  algorithm  obtainable  using 
either  technique  alone,  at  least  for  a  suitable  range  of  n,  m,  U,  and  C.  A  first  result  along  these 
lines  was  obtained  by  Gabow  and  Taijan  [8],  who  developed  an  0(run  log/i  log U  log(nC))-time 
algorithm.  Although  this  time  bound  is  never  less  than  that  of  Goldberg  and  Taijan  [11],  the 
algorithm  does  not  require  sophisticated  data  structures,  whereas  the  Goldberg-Taijan  algorithm 
uses  both  dynamic  trees  [19,20,22]  and  finger  trees  [14,23]. 

Our  work  is  a  continuation  of  efforts  in  this  direction.  We  obtain  an 
0(nm  log  logt/  log  (nC))-t ime  algorithm  for  the  minimum-cost  circulation  problem.  Our  result 
combines  four  known  ideas: 

(1)  Elimination  of  arc  capacities  by  transforming  the  minimum-cost  circulation  problem  into  a 
transportation  problem. 

(2)  Cost  scaling  within  the  e-optimality  framework  as  proposed  by  Goldberg  and  Taijan. 

(3)  A  variant  of  the  Edmonds-Kaip  approach  relying  on  excess  scaling,  as  developed  by  Oriin. 

(4)  The  dynamic  tree  data  structure  of  Sleator  and  Tarjan. 

A  simpler  version  of  our  algorithm  that  does  not  use  the  dynamic  tree  data  structure  runs  in 
0(nm  log U  (1  +  log(nC)/log  logt/))  time.  We  obtain  slightly  better  bounds  for  nonsparse  graphs 
with  very  large  arc  capacities.  We  also  obtain  improved  bounds  for  the  (uncapacitated)  transpor¬ 
tation  problem. 

The  elimination  of  arc  capacities  is  crucial  to  the  efficiency  of  our  algorithms.  An  alterna¬ 
tive  approach  is  to  bound  the  arc  capacities  by  adding  an  extra  outer  capacity-scaling  loop,  as 
suggested  by  Gabow  and  Taijan  [8].  Our  explorations  of  this  approach  lead  to  algorithms  with 
time  bounds  worse  than  those  mentioned  above,  but  the  analytical  methods  we  develop  are  of 
independent  interest  Our  results  using  this  approach  are  described  toward  the  end  of  die  paper. 
For  example,  we  obtain  a  polynomial-time  algorithm  for  the  problem  that  uses  a  classical  net¬ 
work  simplex  algorithm  inside  a  scaling  loop. 

This  paper  consists  of  six  sections  in  addition  to  the  introduction.  In  Section  2  we  define 
the  transportation  problem  and  discuss  its  relationship  with  the  minimum-cost  circulation  prob¬ 
lem.  In  Section  3  we  develop  a  generic  algorithm  for  the  transportation  problem  based  cm  cost 
scaling  and  e-optimality.  In  Section  4  we  refine  the  generic  algorithm  to  use  excess  scaling,  and 
we  analyze  the  resulting  method.  In  Section  5  we  add  the  use  of  dynamic  trees.  In  Section  6  we 
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consider  the  use  of  capacity  bounding  as  an  alternative  way  of  dealing  with  arc  capacities.  In 
Section  7  we  summarize  our  results,  comment  on  the  possible  practicality  of  our  algorithms,  and 
mention  some  open  problems. 

2.  The  Transportation  Problem 

The  minimum-cost  circulation  algorithms  we  develop  in  Sections  3-5  will  be  stated  in  terms 
of  a  related  problem,  the  transportation  problem.  In  order  to  discuss  this  problem,  we  need  some 
terminology.  We  call  C  bipartite  if  V  can  be  partitioned  into  two  sets  5  and  T 
(J  u  T  =  V,  S  r>  T  =  b)  such  that  every  arc  has  exactly  one  vertex  in  S  and  one  in  T.  We  call  ver¬ 
tices  in  S  sources  and  those  in  T  sinks',  we  denote  by  n  i  and  the  sizes  of  S  and  T,  respectively. 
We  call  a  bipartite  network  uncapacitated  if  u(v,w)  *  »  for  each  arc  (v,w)  with  v  e  S  and 
«(v, w)  =  0  for  each  arc  (v,w)  with  v  e  T.  A  supply-demand  vector  d  on  a  bipartite  network  is  a 
mapping  from  V  to  the  real  numbers  such  that  d(v)  £  0  if  v  e  5  (-d(v)  is  the  supply  at  vertex  v), 

d(v)  £  0  if  v  e  T  (d(v)  is  the  demand  at  vertex  v),  and  £  d(v)  =  0  (total  supply  equals  total 

vc  v 

demand).  Given  an  uncapacitated  bipartite  network  C  and  a  supply-demand  vector  d,  the  tran¬ 
sportation  problem  is  that  of  finding  a  minimum-cost  pseudoflow/ satisfying  the  following  con¬ 
straint: 

bj  (v)  =  d(v)  for  all  v  e  V  (supply-demand  constraint)  (6) 


We  call  a  pseudoflow  feasible  if  it  satisfies  (6).  We  call  a  transportation  problem  feasible  if 
it  has  some  feasible  pseudoflow.  Checking  the  feasibility  of  a  traiisportation  problem  can  be  done 
using  any  maximum  flow  algorithm,  e.g.  II  ,2,9,10]. 

There  is  a  well-known,  simple  transformation  that  will  convert  any  minimum-cost  circula¬ 
tion  problem  into  an  equivalent  transportation  problem  [18,24).  Given  a  network  C  =  (V,£)  we 
construct  another  network  G'«(Vu  E,A),  where  A  contains  arcs  ((v,w),v),  ((v,w),w)  and  their 
reversals  for  every  arc  (v.w)  e  E.  The  arcs  ((v,w),v)  and  ((v,w),w)  have  infinite  capacity:  their 
reversals  have  zero  capacity.  Arc  ((v,w),v)  has  cost  zero  and  arc  ((v.wO.w)  has  cost  c(v,w).  We 
define  a  supply-demand  vector  don  G'  by 


d((v,w))  =  -k(v, w)  for  all  (v.w)  c  E,  (7) 

d(v)  -  £  u(v,w)  for  all  v  e  V. 

m:(r.w)t  E 
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Any  circulation  f  on  G  corresponds  to  a  feasible  pseudoflow  f  on  G'  such  that  costff)  = 
cost(f),  given  by  f  ((v,w),w)=/(v,w),  f  ((v,w),v)  =  u(y, w)-/(v.w)  for  each  arc  (v,w)e  E. 
This  correspondence  is  invertible.  Thus  a  solution  to  the  transportation  problem  on  C  gives  a 
solution  to  the  minimum-cost  circulation  problem  on  G.  Observe  that  if  we  regard  £  as  being  the 
set  of  sources  of  G'  and  V  as  being  the  set  of  sinks,  C'  has  n\  -m,  and  n2  =  n;  arc  set  A  has  size 
Am.  (Each  source  has  two  incoming  and  two  outgoing  arcs.) 

We  shall  derive  time  bounds  for  the  transportation  problem  and  translate  them  into  time 
bounds  for  the  minimum-cost  circulation  problem  based  on  the  above  transformation. 

3.  A  Generic  Algorithm  for  the  Transportation  Problem 

We  obtain  a  generic  algorithm  for  the  transportation  problem  by  translating  the  minimum- 
cost  circulation  algorithm  of  Goldberg  and  Tarjan  [9,11]  into  the  setting  of  the  transportation 
problem.  We  modify  the  algorithm  to  be  an  augmenting  path  method;  the  time  bounds  we  derive 
depend  on  this  modification.  We  omit  proofs  of  many  of  the  basic  results,  since  they  are  direct 
translations  of  the  proofs  of  Goldberg  and  Tarjan. 

Let  G  =  (V  =  S  u  7",  A)  be  an  uncapacitated  bipartite  network  with  source  set  S  of  size  n  { , 
sink  set  T  of  size  n2.  arc  set  A  of  size  m.  and  supply-demand  vector  d.  We  denote  the  total  size  of 
V,  i.e.,  ft]  •+  n2,  by  n,  and  min  (/ti,n2]  by  n0.  We  denote  by  V  the  maximum  supply,  i.e.,  U  = 
max  [-d(v)  I  v  e  S) ,  and  by  C  the  maximum  absolute  value  of  an  arc  cost.  Note  that  U  and  C  are 
defined  so  that  the  transformation  of  Section  2  from  a  minimum-cost  circulation  problem  to  a 
transportation  problem  preserves  the  values  of  U  and  C. 

For  a  pseudoflow  / on  G,  we  define  the  excess  e/(v)  of  a  vertex  v  by 


</(v)  =  bj(v)  -d(v). 


(8) 


Thus  a  pseudoflow  is  feasible  if  every  vertex  has  zero  excess.  We  shall  assume  that  the  transpor¬ 
tation  problem  to  be  solved  is  feasible,  i.e.,  there  is  some  feasible  pseudoflow. 

Given  a  pseudoflow  /,  the  residual  capacity  of  an  arc  (v,w)  with  respect  to  /is 


m/v,w)  =  «(v,w) 


(9) 


Arc  (v,w)  is  unsaturated  if  u/y,w)  >  0  and  saturated  otherwise. 

A  price  function  p  on  G  is  a  real-valued  function  on  the  vertices.  Given  a  price  function  p, 
the  reduced  cost  of  an  arc  (v.w)  is 
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Cp(v.w)  =  c(v,w)  +  p(v) -p(w).  (10) 

Let  e  >  0,  let  /  be  a  pseudoflow,  and  let  p  be  a  price  function.  Pseudoflow  /  is  e-optima] 
with  respect  to  price  function  p  if 


cp  (v.w)  £  -e  for  every  unsaturated  arc  (v.w)  (e-opu'mality  constraint).  (1 1) 


Pseudoflow  /is  optimal  with  respect  to  p  if  it  is  e -optima]  for  e  =  0.  The  following  theorem  is  a 
classical  result  of  network  flow  theory  and  follows  from  the  duality  theorem  of  linear  program¬ 
ming. 

Theorem  3.1  [11].  A  feasible  pseudoflow  is  of  minimum  cost  if  and  only  if  it  is  optima]  with 
respect  to  some  price  function  p. 

As  Bertsekas  [4]  discovered,  a  weaker  condition  suffices  if  all  arc  costs  are  integers: 

Theorem  32.  If  all  arc  costs  are  integers  and  E  <  then  a  feasible  flow  is  of  minimum  cost  if 

2  n0 

and  only  if  it  is  e-optimal  with  respect  to  some  price  function  p. 

Proof.  Analogous  to  the  proof  or  Theorem  2.3  of  [11],  using  the  fact  that  G  is  bipartite  and  hence 
any  simple  cycle  contains  at  most  2 n0  vertices.  □ 

In  the  remainder  of  this  paper  (except  in  some  of  the  concluding  remarks  of  Section  7),  we 
shall  assume  that  all  arc  costs  are  integers;  thus  Theorem  3 2  applies. 

Our  algorithm  applies  cost  scaling  based  on  Theorem  32.  It  uses  a  cost-scaling  factor 
ki  2.  It  maintains  a  price  function  p  and  an  error  parameter  e.  Initially  e  =  C  and  p  is  identi¬ 
cally  zero.  The  algorithm  consists  of  repeating  the  following  step  until  the  termination  condition 
is  satisfied. 

Cost-Scaling  Step.  Let  fix  the  identically  zero  pseudoflow.  By  modifying  f  and  p,  find  a  feasi¬ 
ble  pseudoflow  f  and  a  price  function  p'  such  that  f  is  e-optimal  with  respect  to  p'.  If 

c  <  — — ,  stop.  Otherwise,  let  p  be  defined  by  p(v)  =  p'(v)  +  e  if  v  e  S,  p(v)  =  p'(v)  if  v  e  T; 
2/Jo 

replace  e  by  tik. 
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Note  that  only  the  price  function  p  is  carried  over  from  iteration  to  iteration;  the  pseudoflow 
is  ret  to  zero  after  each  iteration. 

Lemma  33.  At  the  beginning  of  a  cost-scaling  step,  /  (the  zero  pseudoflow)  is  e-optimum  with 
respect  to  p. 

Proof.  Any  arc  (v,w)  that  is  unsaturated  with  respect  to  /has  v  e  S  and  w  e  T.  Suppose  that  the 
current  cost-scaling  step  is  the  first  Then  cp(v,w)  =  c(v,w)  £  -C  «  -£-  Suppose  on  the  other 
hand  that  the  cost  scaling  step  is  not  the  first  Let  f  be  the  pseudoflow  and  p'  the  price  function 
computed  in  the  previous  step.  Then  (v.w)  is  unsaturated  with  respect  to  f ,  since  u(v,w)  = «», 
but  any  pseudoflow  has  all  arc  flows  finite.  Thus  cp(y,w)  =  cy  (v,w)  +  kt  £  -kt  ♦  ke  2  0,  since  f 
is  kt  optima]  with  respecttop',andp(v)  =  />'(v)  +  *£P(*')=p,(H')-  D 

Theorem  3.4.  The  transportation  algorithm  is  correct  and  terminates  after  O  (1  +  log*(n0  C)) 
iterations. 

Proof.  Correctness  follows  from  Theorem  3.2.  The  bound  on  the  number  of  iterations  is  obvi¬ 
ous.  □ 

The  heart  of  the  algorithm  is  the  conversion  of  an  e -optimal  pseudoflow  into  an  e -optimal 
feasible  pseudoflow  and  the  corresponding  modification  of  the  price  function.  We  call  this  the 
refinement  computation.  Our  generic  refinement  algorithm  consists  of  a  sequence  of  two  kinds  of 
local  transformations,  one  of  which  modifies  the  pseudoflow  and  the  other  of  which  modifies  the 
price  function.  To  define  the  transformations,  we  use  the  following  terminology.  A  vertex  v  is 
active  if  ej  (v)  >  0.  An  unsaturated  arc  (v,w)  is  eligible  if  cp  (v,w)  <  0.  The  refinement  algorithm 
consists  of  repeating  the  following  steps,  in  any  order,  until  no  vertex  is  active,  and  then  defining 
/*/.  p'*=p: 

push  (v.w); 

Applicability:  Vertex  v  is  active,  u/  (v.w)  >  0,  and  c,  (v.w)  S  e. 

Action:  Push  up  to  8  =  min{e/  (v).  uj  (v.w) )  units  of  flow  from  v  to  w  by  increasing/  (v.w) 
by  an  amount  up  to  6. 


relabel  (v). 

Applicability:  Vertex  v  is  reachable  from  some  active  vertex  by  a  path  of  eligible  arcs,  and 
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there  is  no  eligible  arc  (v.h'). 

Action:  Replace p(v )  by  max  {  p(w)  -c(v, w)  -e). 

Lemma  35.  Any  pushing  or  relabeling  step  preserves  e-optimality.  A  relabeling  of  a  vertex  v 
decreases  p(v)  by  at  least  e. 

Proof.  Analogous  to  the  proofs  of  Lemmas  5.2  and  5.3  of  [1 1].  □ 

Lemma  3.6.  The  price  of  a  vertex  v  decreases  by  0  (knQt)  during  refinement.  Hence  v  is  rela¬ 
beled  O(kn0)  times. 

Proof.  Analogous  to  the  proof  of  Lemmas  5.7,  5.8,  and  5.9  of  [11].  The  bound  on  price  changes 
in  Lemma  5.7  of  [1 1]  is  O  (ne),  where  the  cost  scaling  is  by  a  factor  of  two.  Revising  the  argu¬ 
ment  to  include  a  cost  scaling  factor  of  k  yields  an  0  (Lue)  bound.  Observing  that  G  is  bipartite, 
and  hence  that  any  simple  path  in  G  contains  at  most  2n0  arcs,  reduces  the  bound  to  O  ( kn0t ).  □ 

Now  we  describe  a  version  of  the  refinement  algorithm  that  is  based  on  the  idea  of  finding 
augmenting  paths.  The  algorithm  uses  a  fixed  incidence  list  /(v)  for  each  vertex  v.  This  list  con¬ 
tains  each  arc  (v.w).  One  such  arc  is  designated  the  current  arc  out  of  v.  Initially  the  current  arc 
out  of  v  is  the  first  arc  on  /(v).  The  algorithm  repeatedly  attempts  to  find  a  path  of  eligible  arcs 
from  an  active  vertex  to  a  vertex  of  negative  excess.  When  such  a  path  is  found,  flow  is  pushed 
along  it.  To  find  such  paths,  the  algorithm  uses  depth-first  search,  implemented  using  a  stack  S. 
During  a  search,  vertices  are  relabeled  as  necessary  to  extend  the  path.  The  algorithm  consists  of 
initializing  5  to  be  empty  and  repeating  the  following  steps  until  termination  occurs  in  Step  1. 

Step  1  (stan  new  path).  If  there  are  no  active  vertices,  stop.  Otherwise,  select  some  active  venex 
v  and  push  it  onto  5.  Co  to  Step  2. 

Step  2  (extend  path).  Let  v  be  the  top  vertex  on  S.  While  the  current  arc  of  v  is  not  eligible, 
replace  the  current  arc  by  the  next  arc  on  /(v).  If  the  end  of  /(v)  is  reached  without  finding  an  eli¬ 
gible  are,  go  to  Step  3.  If  an  eligible  are  (v,w)  is  found,  test  whether  e/(w)  £  0.  If  so,  push  w 
onto  S  and  repeat  Step  2;  if  not,  go  to  Step  4. 

Step  3  (relabel).  Relabel  v,  the  top  vertex  on  5.  Reset  the  current  arc  of  v  to  be  the  first  arc  on 
/(v).  If  v  is  not  the  only  vertex  on  S,  pop  it  from  S.  Go  to  Step  2. 
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Step  4  (augment).  Let  6  be  any  positive  quantity  not  more  than  the  minimum  of  e/  (v) 
and  min  { u/(x,y )  I  x  is  on  S,  and  (x,y)  is  the  current  arc  out  of  x }.  For  each  current  arc  (x,y)  such 
that  x  is  on  S ,  increase  f(x,y)  by  6.  Empty  S  and  go  to  Step  1. 

We  call  this  method  the  augmenting  path  version  of  the  refinement  algorithm,  or  the  aug¬ 
menting  path  algorithm  for  short.  We  call  an  execution  of  Step  4  an  augmentation. 

Lemma  3.7.  The  augmenting  path  algorithm  maintains  the  invariant  that  there  is  no  cycle  of  eli¬ 
gible  arcs. 

Proof.  Analogous  to  the  proof  of  Corollary  5.6  of  [1 1  ].  □ 

Remark..  The  proof  of  Lemma  3.7  uses  the  fact  that  pushes  take  place  only  along  eligible  arcs.  □ 
Lemma  3.8.  The  maximum  size  of  S  is  at  most  2no ■ 

Proof.  The  vertices  on  S  always  define  a  path  of  eligible  arcs.  By  Lemma  3.7  such  a  path  is  sim¬ 
ple.  The  fact  that  G  is  bipartite  gives  the  claimed  bound  on  the  size  of  S.  □ 

Theorem  3.9.  The  augmenting  path  algorithm  is  correct  and  runs  in  O(kn0m)  time  plus  0(n 0) 
time  per  augmentation. 

Proof.  Correctness  follows  from  Lemma  3.5.  We  bound  the  running  time  as  follows.  The 
number  of  additions  to  S  equals  the  number  of  pops  from  5.  The  number  of  pops  from  5  is  O(n0 ) 
per  augmentation  by  Lemma  3.8  plus  at  most  one  per  relabeling.  The  time  to  relabel  a  vertex  v  is 
0(  I  /(v)  I),  which  is  also  the  time  spent  in  Step  2  changing  current  arcs  of  v  between  relabelings. 
By  Lemma  3.6,  the  relabeling  time  and  time  spent  changing  current  arcs,  summed  over  all  ver¬ 
tices,  is  O(kn0m).  An  execution  of  Step  2  that  does  not  change  the  current  arc  of  v  causes  an 
addition  to  S.  The  time  to  do  an  augmentation  is  O(no).  The  claimed  time  bound  follows.  □ 

4.  Bounds  for  the  Augmenting  Path  Algorithm 

In  this  section  we  derive  time  bounds  foT  various  versions  of  the  augmenting  path  algo¬ 
rithm.  Observe  that  there  are  two  kinds  of  freedom  in  this  algorithm,  in  the  choice  of  starting  ver¬ 
tices  for  augmenting  paths  in  Step  1,  and  in  the  amount  by  which  the  flow  is  augmented  in  Step  4. 

Let  us  first  analyze  the  simple  method  in  which  each  augmentation  is  by  an  amount  that  is 
as  large  as  possible;  that  is,  in  Step  4,  6  is  selected  as  follows:  6  =  min  {  ej(v),  min  (u/x.y)  I  x  is 
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on  5,  x  *  *,  and  (x.y)  is  a  current  arc} ).  With  this  method,  each  augmentation  either  reduces  the 
number  of  active  vertices  by  one,  reduces  the  number  of  vertices  of  negative  excess  by  one,  or 
saturates  an  arc.  Lemma  3.6  implies  that  the  total  number  of  arc  saturations  is  O(kn0m)  (see 
Lemma  5.10  of  [1 1]);  hence  the  total  running  time  of  the  augmenting  path  algorithm  is  O(knlm), 
or  O(nlm)  if  k  is  chosen  equal  to  two.  This  bound  is  analogous  to  Dinic's  bound  of  O  ( n2m )  for 
the  maximum  flow  problem  [5];  indeed,  the  augmenting  path  algorithm  itself  can  be  viewed  as  an 
analogue  of  Ahuja  and  Orlin’s  variant  [17j  of  Dime's  algorithm. 

We  obtain  a  better  bound  (if  all  arc  capacities  are  integral  and  not  too  large)  by  using  excess 
scaling.  This  method  is  based  on  Orlin’s  variant  116]  of  the  capacity-scaling  algorithm  of 
Edmonds  and  Karp  [6]  for  the  minimum-cost  circulation  problem  and  is  also  analogous  to  the 
maximum  flow  algorithm  of  Gabow  [7],  Henceforth  (except  in  Section  7)  we  shall  assume  that 
all  arc  capacities  are  integral. 

The  excess-scaling  algorithm  maintains  an  estimate  A  of  the  maximum  excess.  Initially  A 
is  the  largest  power  of  two  not  exceeding  U.  The  algorithm  maintains  two  invariants: 

(i)  The  sum  of  all  positive  excesses  is  at  most  2nA; 

(ii)  The  residual  capacity  of  any  arc  is  either  infinity  or  an  integer  (possibly  zero)  multiple  of 
A. 


In  Step  1 ,  the  algorithm  always  chooses  a  starting  vertex  v  with  tj  (v)  £  A;  if  no  such  vertex 
exists,  the  algorithm  replaces  A  by  A/2  and  tries  again.  In  Step  4,  the  algorithm  always  pushes  A 
units  of  flow  along  the  augmenting  path.  The  choice  of  starting  vertices  guarantees  that  invariant 
(i)  is  maintained;  immediately  after  A  changes,  the  sum  of  all  positive  excesses  is  at  most  2nA. 
Augmenting  by  A  preserves  invariant  (ii),  which  in  turn  guarantees  that  A  units  of  flow  can  actu¬ 
ally  be  pushed  each  time  an  augmentation  occurs.  When  A  =  1/2,  all  excesses  are  zero,  and  the 
algorithm  terminates.  We  call  a  maximal  period  of  time  during  which  A  stays  constant  a  phase  of 
the  algorithm. 

Lemma  4.1  The  total  number  of  augmentations  done  by  the  excess-scaling  algorithm  is 
0(n  log  U). 

Proof.  Each  augmentation  either  reduces  the  number  of  vertices  with  negative  excess  by  one  or 
reduces  the  sum  of  positive  excesses  by  A.  By  (i),  the  latter  case  can  occur  only  O(n)  times  dur¬ 
ing  a  phase  The  number  of  phases  is  OGog  U).  The  bound  follows.  □ 
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Theorem  42.  The  excess-scaling  version  of  the  augmenting  path  algorithm  runs  in 
O(n0(km  +  nlog  U))  time.  Using  this  method  in  the  transportation  algorithm  gives  a  bound  for 
the  transportation  problem  of  0{nQ{km  +  n  log  U)  logt  (nC))  time  for  any  k  such  that  2  S  k  £  nC. 

Choosing  k  =  min  {2+  —  log  U,  nC  }  yields  the  following  time  bounds  for  the  transportation 

m 

problem: 


O(nom  log  ( nC ))  if  log  U  <  2m!n\ 


O(n0n  log U  (1  +log  (nC)/log  (—  log  U)))  if  log  U  2  2 min. 

m 


Proof.  Immediate  from  Lemma  4.1  and  Theorem  3.4.  □ 

Corollary  43.  The  excess-scaling  version  of  the  transportation  algorithm  combined  with  the 
transformation  of  Section  2  will  solve  a  minimum-cost  circulation  problem  in 
0(nm  log  U  (1  +  log  (nC)/ log  log  U ))  time. 

By  changing  the  excess-scaling  algorithm  slightly,  wc  can  obtain  a  bound  of 
0(n  o  log  U  +  n  log  min  {  n,  U ) )  on  the  number  of  augmentations.  This  is  an  improvement  on  the 
bound  of  Lemma  4.1  only  if  log  U  =  to((n/no)logn)>  which  only  holds  if  U  grows  nonpolynomi- 
ally  with  n.  Nevertheless,  we  shall  present  the  result,  since  it  suggests  the  possibility  of  obtaining 
an  O  (n  log n)  bound  on  the  number  of  augmentations  for  some  suitable  modification  of  the  algo¬ 
rithm.  We  shall  assume  that  n  1 2t  n2,  >  c.  n0  =  n2,  which  is  without  loss  of  generality:  if  n,  <  n2, 
exchange  the  source  set  and  the  sink  set  and  negate  the  supply-demand  vector. 

We  modify  the  excess-scaling  algorithm  by  changing  Step  I  to  the  following: 

Step  I If  there  are  no  active  vertices,  stop.  Otherwise,  if  some  active  vertex  v  c  5  has  an  outgo¬ 
ing  arc  (v.w)  such  that  /  (v,w)  >  7n&,  increase  the  flow  on  (v,w)  by  e/(v)  and  repeat  Step  1 
(We  call  this  a  special  push.)  Otherwise,  if  some  active  vertex  v  has  ey  (v)  2  A,  push  vertex  v 
onto  S  and  go  to  Step  2.  Otherwise,  replace  A  by  A/2  and  repeal  Step  1 '. 

Before  analyzing  the  modified  algorithm  in  detail,  we  make  several  observations.  A  special 
push  is  actually  a  push,  since  if  /(v.w)  >  0  then  uy  (w,v)  >  0,  which  implies  by  e-optimality  that 
cp  (w,  v)  2  -e  and  by  cost  antisymmetry  that  cp  (v.w)  S  e.  A  special  push  maintains  the  invariant 


-12- 


thai  there  is  no  cycle  of  eligible  arcs.  Once  a  vertex  v  e  S  has  zero  excess,  its  excess  remains 
zero  until  the  end  of  the  algorithm.  The  excess  on  any  vertex  v  e  S  is  strictly  less  than  2A.  The 
total  flow  moved  by  special  pushes  is  thus  less  than  2n  A.  The  total  flow  moved  during  augmen¬ 
tation:  that  decrease  the  number  of  vertices  with  negative  excess  is  at  most  n  A.  The  total  flow 
moved  by  other  augmentations  during  a  single  phase  is  at  most  the  sum  of  the  positive  excesses, 
which  is  at  most  2nA.  Thus  the  total  flow  moved  from  a  given  time  until  the  end  of  the  algorithm 

m 

is  at  most  3/iA  +  J  2/»A/2‘  =  7nA.  It  follows  that  once  an  arc  (v,w)  has  flow  exceeding  7nA,  its 

i-0 

flow  remains  positive  until  the  end  of  the  algorithm,  and  (w,v)  can  never  be  saturated.  We  call  an 
arc  (v,w)  that  can  never  be  saturated  open.  The  modified  algorithm  maintains  invariant  (i)  (the 
sum  of  all  positive  excesses  is  at  most  2nA)  and,  in  place  of  invariant  (ii),  the  following: 

(ii)  ’  Every  arc  (v,w)  is  either  open  or  has  a  residual  capacity  that  is  an  integer  ( possibly  zero) 
multiple  of  A. 

We  can  verify  invariant  (ii) '  by  induction  on  the  number  of  steps  taken  by  the  algorithm, 
simultaneously  showing  that  every  augmentation  in  Step  4  can  actually  move  A  units  of  flow. 

Lemma  4.4.  The  total  number  of  augmentations  made  by  the  modified  excess-scaling  algorithm 
isO  (n0logl!  +  n  logmin(  n,i/}). 

Proof.  Consider  a  vertex  v  e  S.  There  is  at  most  one  augmentation  starting  from  v  per  phase. 
Suppose  that  the  first  augmentation  from  v  is  during  phase  i.  This  augmentation  moves  A  units  of 
flow.  Henceforth  until  the  end  of  die  algorithm  there  is  always  an  arc  (v.w)  with  /  (v.w)  2  A  in. 
After  2  log  n  +  3  more  phases,  the  current  value  of  the  excess  estimate  is  A'  =  A/8n2,  and  there  is 
some  arc  (v.w)  with  /(v.w)  2t  8n  A'.  When  such  an  arc  exists,  if  not  before,  the  excess  at  v  is 
reduced  to  zero  by  a  special  push.  Hence  v  can  have  positive  excess  only  during 
O  (logmin{  n,U))  phases,  and  there  are  0(n  log  min(  n,l/))  augmentations  starting  from  ver¬ 
tices  in  5. 

Now  consider  a  vertex  v  €  T.  If  v  does  not  receive  additional  flow  from  special  pushes  dur¬ 
ing  a  phase,  there  can  be  at  most  one  augmentation  starting  from  v  during  the  phase.  A  special 
push  to  v  moves  less  than  2A  units  of  flow  to  v.  which  can  account  for  at  most  two  augmentations 
starting  from  v  during  the  phase.  We  charge  such  augmentations  to  die  corresponding  special 
push.  Since  there  are  only  n  special  pushes,  the  number  of  augmentations  starting  from  vertices 
in  T,  summed  over  all  phases,  is  O  (n  +  n o  log  U ).  This  gives  the  desired  bound.  □ 
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In  presenting  time  bounds  for  the  modified  method,  we  assume  that 
log  V  =  fl((n  /n0)  log  n),  since  otherwise  the  bounds  are  the  same  as  those  in  Theorem  4.2  and 
Corollary  4.3. 

Theorem  4.5.  Assume  that  log  U  -  fi((n/n0)logn).  Then  the  modified  excess-scaling  version 
of  the  augmenting  path  algorithm  nms  in  O(n0  (bn  +  n0  log  U))  time.  Using  this  method  in 
the  transportation  algorithm  gives  a  time  bound  for  the  transportation  problem  of 
O(n0(bn  +  n0  log  U)  log*  (n0C))  for  any  k  such  that  2  £k& n0C.  The  choice  of 

n0 

k  =  min  {  2  +  —  log  U,  noC)  yields  the  following  time  bounds  for  the  transportation  problem: 


O(nom  log  («oC))  if  log  U  <  2m  / no’. 


0(n o  log  V  (1  +  log  (n0C)l  log  ( —  log  £/)))  if  log  U  2  2m/n0. 

m 


Corollary  4.6.  If  log  U  =  ft  (—  logn),  the  modified  excess-scaling  version  of  the  transportation 

n 

algorithm  combined  with  the  transformation  of  Section  2  will  solve  a  minimum-cost  circulation 
problem  in  O  (n2  log  U(1  +  log  ( nC )/  log  log U))  time. 

Remark.  Every  bound  derived  in  this  Section  remains  valid  if  each  occurrence  of  the  parameter  V 

is  replaced  by  another  smaller  parameter  U* .  For  the  transportation  problem, 

U*=4+  (~d(v))/n.  For  the  minimum-cost  circulation  problem,  17*  =4+  £  u(v,w)im. 

*«  S  (»,»■)«  E 

The  bound  on  the  number  of  augmentations  in  Lemma  4.1  can  be  reduced  to  O  (n  log V*)  by 

observing  that  the  sum  of  positive  excesses  is  initially  at  most  nU *,  which  implies  that  the 

number  of  augmentations  dunng  phases  in  which  A  >  U*  is  O  (n).  Similarly,  die  bound  in 

Lemma  4.4  can  be  reduced  to  O  (no  log U*  +n  log  min  Corresponding  improvements 

in  the  bounds  of  Theorem  4.2,  Corollary  4.3,  Theorem  4.5,  and  Corollary  4.6  follow.  These 

improved  bounds  are  analogous  to  the  bound  Edmonds  and  Karp  obtained  for  their  transportation 

algorithm  [6]. 
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5.  Use  of  Dynamic  Trees 

The  algorithms  discussed  in  Section  4  ait  quite  simple  and  do  not  require  the  use  of  any 
complicated  data  structures.  By  adding  the  use  of  dynamic  tree*'  [19,20,22],  we  can  improve  the 
bounds  derived  in  Section  4  by  almost  a  logarithmic  factor.  Our  use  of  dynamic  trees  is  analo¬ 
gous  to  their  use  in  other  network  flow  algorithms  [2,9,10,11,12,19,22]. 

The  dynamic  tree  data  structure  allows  the  maintenance  of  a  collection  of  vertex-disjoint 
rooted  trees,  each  arc  of  which  has  an  associated  value.  Each  tree  is  an  in-tree;  that  is,  if  vertex  v 
is  a  child  of  vertex  w,  there  is  a  tree  arc  from  v  to  w.  Each  vertex  in  a  tree  is  regarded  as  being 
both  an  ancestor  and  a  descendant  of  itself.  The  data  structure  supports  the  following  seven 
operations: 

find-roo^v):  Find  and  return  the  root  of  the  tree  containing  vertex  v. 

fwd-size(vy.  Find  and  return  the  number  of  vertices  in  the  tree  containing  vertex  v. 

find-value(v ):  Find  and  return  the  value  of  the  tree  arc  leaving  v.  If  v  is  a  tree  root,  the 
value  returned  is  infinity. 

find-min(y):  Find  and  return  the  ancestor  w  of  v  with  find-valiuKyv)  minimum.  In  case 
of  a  tie,  choose  the  vertex  h>  closest  to  the  tree  root. 

change-value(v,x):  Add  real  number  x  to  the  value  of  every  arc  along  the  path  from  v 
to  find-root(v). 

Unk(y,w,x)'.  Combine  the  trees  containing  v  and  w  by  making  w  the  parent  of  v  and  giving 
the  arc  (v,w)  the  value  x.  This  operation  does  nothing  if  v  and  w  are  in  the  same  tree  or 
if  v  is  not  a  tree  root. 

cu/(v):  Break  the  tree  containing  v  into  two  trees  by  deleting  the  edge  joining  v  to  its 

parent;  return  the  value  of  the  deleted  edge.  This  operation  breaks  no  edge  and  returns 
infinity  if  v  is  a  tree  root 

A  sequence  of  /  tree  operations,  starting  with  an  initial  collection  of  singleton  trees,  takes 
0(1  log  (r  +  1))  time  if  z  is  the  maximum  tree  size  [J0,19,20,22J. 

We  use  this  data  structure  to  represeni  a  subset  of  the  eligible  current  arcs.  The  value  of  an 
arc  is  its  residual  capacity.  The  data  structure  allows  flow  to  be  moved  along  an  entire  path  ai 
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once,  rather  than  along  one  arc  at  a  time. 

In  applying  this  data  structure  to  the  transportation  problem,  we  can  improve  the  resulting 
time  bounds  if  we  take  advantage  of  the  special  structure  of  the  problem,  specifically  the  fact  that 
G  is  bipartite,  and  hence  so  is  every  dynamic  tree.  Let  us  assume  that  n  j  i  f»2,  i.e.  I  S  !  £  I T  I. 
We  redefine  the  size  of  a  dynamic  tree  to  be  the  number  of  vertices  of  T  it  contains.  This  changes 
the  semantics  of  the  find  size  operation,  but  does  not  affect  its  implementation  significantly.  We 
also  modify  the  data  structure  so  that  any  dynamic  tree  contains  at  most  twice  as  many  vertices  as 
its  size.  To  do  this  we  introduce  an  extra  layer  of  abstraction.  We  represent  each  of  the  actual 
dynamic  trees  (the  ones  manipulated  by  the  operations)  by  a  virtual  dynamic  tree,  which  consists 
of  the  actual  tree  with  all  leaves  in  5  deleted.  Each  of  the  deleted  leaves  has  a  pointer  to  its  parent 
in  the  actual  tree,  and  has  stored  with  it  the  value  of  the  outgoing  tree  arc.  Every  virtual  tree  con¬ 
tains  a  number  of  vertices  at  most  twice  its  size,  since  every  virtual  tree  vertex  in  S  has  a  virtual 
tree  child  in  T.  Every  operation  on  actual  trees  translates  into  0(1)  operations  on  virtual  trees.  It 
follows  that  a  sequence  of  /  operations  on  actual  dynamic  trees  takes  0(1  log  (z  +  1))  time,  where 
2  is  the  maximum  tree  size  according  to  the  new  definition  of  size. 

The  following  version  of  the  excess-scaling  algorithm  uses  these  modified  dynamic  trees. 
In  addition  to  an  excess  estimate  A,  the  algorithm  uses  a  fixed  bound  z,  15zSn0.onlhe  max¬ 
imum  size  of  a  dynamic  tree.  The  algorithm  maintains  a  stack  S  that  defines  a  path  of  eligible 
current  arcs  as  follows:  if  vertex  v  appears  just  below  vertex  w  on  S,  then  the  tree  path  from  v  to 
find-root(v )  followed  by  the  arc  (find-root(v), w)  is  a  path  of  eligible  current  arcs.  Initially  S  is 
empty  and  each  vertex  forms  a  one-vertex  dynamic  tree.  The  algorithm  consists  of  repeating  the 
following  steps  until  termination  occurs  in  Step  1. 

Step  1  (start  new  path).  If  no  vertex  has  positive  excess,  stop.  Otherwise,  if  no  vertex  has  excess 
at  least  A,  replace  A  by  A/2  and  repeat  Step  1.  Otherwise,  let  v  be  a  vertex  of  excess  at  least  A. 
Push  v  onto  S  and  go  to  Step  2. 

Step  2  (extend  path).  Let  v  be  the  top  vertex  on  S.  Compute  w  =  find  root(v).  If  tj  (w)  <  0,  go  to 
Step  4.  Otherwise,  while  the  current  arc  of  w  is  not  eligible,  replace  the  current  arc  of  w  by  the 
next  arc  on  /  (w).  If  the  end  of  /  (w)  is  reached  without  finding  an  eligible  arc,  go  to  Step  3.  If  an 
eligible  arc  (w,x)  is  found,  test  whether  find-size(v )  ♦  find-size(x)  S  z.  If  so,  perform 
link(wa,  u(v,  w)  -/(v,h>)).  If  not,  push  x  onto  S.  Repeat  Step  2. 

Step  3  (relabel).  Relabel  w.  For  each  tree  arc  (y, w),  perform  cut(y).  If  v  =  w  and  v  is  not  the 
only  vertex  on  S,  pop  v  from  S.  Go  to  Step  2. 
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Step  4  (augment).  Add  A  to  e/w). 

Sup  4a.  Perform  change-valueiv,  -A).  While  fuid-value(find-min(y))  -  0,  perform 
cu t  (find-min(v)).  Go  to  Step  4b. 

Step  4b.  Pop  v  from  S.  If  S  is  empty,  subtract  A  from  ej  (v)  and  go  to  Step  1.  Otherwise, 
let  x  =  v  and  replace  v  by  the  new  top  vertex  on  S.  Let  w  =  fmd-root(y).  Add  At of  (h>,x) 
and  go  to  Step 

This  algorithm  stores  flow  explicitly  for  arcs  that  are  not  in  dynamic  trees  and  implicitly  for 
tree  arcs.  Whenever  a  cut  is  performed,  the  arc  cut  must  have  its  flow  restored  to  its  correct 
current  value.  When  the  algorithm  terminates,  every  arc  still  in  a  dynamic  tree  must  have  its 
correct  flow  computed.  These  compulations  have  been  omitted  from  the  description  above. 

The  analysis  of  this  algorithm  is  similar  to  the  analysis  of  other  network  flow  algorithms 
that  use  dynamic  trees,  e.g.  [2,9,10,11,12,19,22].  Since  with  this  method  the  time  bound  for  the 
transportation  problem  is  not  improved  by  using  a  non-constant  cost-scaling  factor,  we  shall 
choose  k  *  2.  The  total  number  of  links  and  cuts  performed  in  the  dynamic  tree  version  of  the 
excess-scaling  algorithm  is  O(n0m)  (see  e.g.  Lemma  7.2  of  [11]),  taking  time 
O(n0m  log  (z  +  1)).  The  proof  of  Lemma  4.1  is  valid  for  this  version  of  the  excess-scaling  algo¬ 
rithm,  which  means  that  there  are  0(n  lo gU)  augmentations.  The  definition  of  the  algorithm 
guarantees  that  if  v  and  w  are  consecutive  vertices  on  S,  and  w  is  not  the  top  vertex  on  S  ,  then 
find-sizeiy)  +  find-size(w)  >  z.  Thus  either  the  tree  containing  v  or  the  tree  containing  w  has  size 
exceeding  z/2.  Since  every  vertex  on  5  is  in  a  different  dynamic  tree,  the  maximum  height  of  S  is 
O(n0lz),  and  the  time  per  augmentation  is  O((n0fz )  log  (z  +  1)).  Thus  we  obtain  the  following 
result 

Theorem  5.1.  The  dynamic  tree  version  of  the  excess -scaling  algorithm  runs  in 
O(n0  (m- f  y  log U)  log  (z  +  1))  time,  for  any  z  satisfying  1  £  r£  «0.  With  this  method,  the  tran¬ 
sportation  algorithm  runs  in  0(«o(w  +  -jlogl/)log(r  +  l)log(noO)  time.  Choosing 
t  *  min  { 1  +  —  lose/,  n  o  )  gives  the  following  time  bounds  for  the  transportation  problem: 


O(n0m  log  (2  +  —  lo  gU)  log  (*oO)  if  log  V  <  n0m/n; 

ffX 
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0(n  \ogU  logno  log  ( noC ))  if  log  U  i  n^mln. 


Remark.  The  bound  in  Theorem  5.1  for  the  case  log  U  £  n0m/n  i*  ot  an  interesting  one,  since  a 
better  bound  of  O(n0m  log  n0  log  (n0O)  can  be  obtained  by  implementing  the  generic  augment¬ 
ing  path  algorithm  (without  excess  scaling)  using  dynamic  trees.  □ 

Corollary  52.  The  dynamic  tree  implementation  of  excess  scaling  can  be  used  to  solve  the 
minimum-cost  circulation  problem  in  0(nm  log  logt/  log  («oC))  time  if  log  V  £  n.  If  log  U  >  n, 
a  bound  of  0(nm  logn  log  ( n0C ))  is  obtainable  with  a  dynamic  tree  implementation  of  the  aug¬ 
menting  path  algorithm  without  excess  scaling. 

We  can  reduce  the  bound  on  the  number  of  augmentations  to 
O(n0  log U  +  n  log  min(n,i/} )  by  modifying  the  excess-scaling  method  as  in  Section  4.  This 
leads  to  the  following  results. 

Theorem  53.  Assume  that  log  U  =  ((«/«o)  logn).  Then  the  dynamic  tree  implementation  of 

n0 

the  modified  excess  scaling  algorithm  runs  in  O(n0  (m  +  —  log  U)  log  (z  +  1))  time,  for  any  z 

satisfying  lSzSno-  With  this  method,  the  transportation  algorithm  runs  in 

n0  «o 

O(n0  (m  +  —  log  U)  log  (z  +  1)  log  (noC))  time.  Choosing  z  =min  {1  +  —  log  U,  n0  J  gives 
z  m 

the  following  lime  bounds  for  the  transportation  problem: 


O(n0m  log  (2  +  —  log  U)  log  ( n0C ))  if  log  U  <  m; 
m 


O(n0  log  U  log  no  log  (noC))  if  log  U2m. 


Corollary  5.4.  If  log  U  *  ft  ( —  logn),  the  minimum-cost  circulation  problem  can  be  solved  in 

H 

0(nm  log  (2  +  —  log  V)  log  ( nC ))  time. 


Remark.  Every  bound  derived  in  this  Section  remains  valid  if  each  occurrence  of  U  is  replaced 
by  £/*,  where  V*  is  as  defined  in  the  remark  at  the  end  of  Section  4.  This  follows  from  die 
corresponding  improvements  in  the  bounds  of  Lemmas  4.1  and  4.4  discussed  in  that  remark.  □ 
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6.  The  Capacity  Bounding  Technique 

The  results  derived  in  Sections  3-5  depend  crucially  on  the  elimination  of  arc  capacities  via 
the  transformation  to  a  transportation  problem  discussed  in  Section  2.  One  may  ask  whether 
there  is  some  more  direct,  or  at  least  alternative,  way  to  deal  with  arc  capacities.  A  question  that 
turns  out  to  be  related  is  whether  the  Ahuja-Orlin  excess-scaling  algorithm  for  the  maximum  flow 
problem  [1,2]  generalizes  in  a  natural  way  to  the  minimum-cost  flow  problem  via  cost  scaling,  in 
analogy  with  the  generalization  of  other  maximum  flow  algorithms  to  this  problem  [9,10). 

In  this  section  we  show  that  the  answer  to  both  of  these  questions  is  a  qualified  “yes.”  We 
consider  the  minimum -cost  circulation  problem  as  defined  in  Section  1,  with  integer  capacities 
and  costs.  We  propose  a  way  of  solving  this  problem  using  an  outer  capacity-scaling  loop  whose 
effect  is  to  convert  die  original  problem  into  a  sequence  of  0(log  U)  problems  in  each  of  which 
the  arc  capacities  are  integers  bounded  by  m.  This  idea  was  used  by  Gabow  and  Taijan  [8].  The 
method  requires  a  standard  “vertex-splitting"  transformation.  To  solve  the  resulting  capacity- 
bounded  problems,  we  propose  a  modification  of  die  Ahuja-Oriin  excess-scaling  maximum  flow 
algorithm  nested  inside  an  e-scaling  loop.  The  resulting  triple  scaling  algorithm  runs  in 
0((n2logm  +  nm)  log U  log (nC))  time  using  no  fancy  data  structures.  Although  this  algorithm 
has  an  inferior  complexity  bound  as  compared  to  the  bounds  obtained  in  previous  sections,  the 
method  and  its  analysis  have  independent  interest 

Now  we  give  details.  The  outer  capacity-scaling  loop  constructs  two  minimum-cost  circu¬ 
lation  problems  at  each  iteration,  a  target  problem  and  a  restricted  ( capacity-bounded )  problem. 
These  are  obtained  from  the  current  network  as  follows.  First  the  next  bit  of  precision  is  added 
»o  the  arc  capacities  by  doubling  the  current  capacity  and  adding  one  to  the  capacity  of  each  arc 
(v,w)  such  that  the  current  bit  of  u(y,w)  is  1.  We  denote  the  resulting  capacity  function  by  u' . 
The  target  problem  is  ( C,u',c ),  where  C  is  the  original  graph  and  c  is  the  original  cost  function. 

The  restricted  problem  is  obtained  from  the  target  problem  by  bounding  the  flow  through 
every  vertex  of  the  target  problem  by  m.  More  formally,  the  restricted  problem  is  obtained  by 
splitting  each  vertex  v  of  C  *  (V,£)  into  two  to  obtain  the  graph  G'  *  (V\£'),  where  V'  contains 
vertices  Vi  and  v2  for  each  v  €  V,  and  £'  contains  an  arc  (vt,v2)  of  capacity  u(vi,v2)  =  m  and 
cost  c(v(,v2)  =  0  for  each  v  e  V  and  an  arc  (v2,wi)  of  capacity  «"(v2,h',)  =  «'(v,w)  and  cost 
c'(v2, W|)  =  c(v,w)  for  each  arc  (v,w)  e  £.  Network  G'  also  contains  opposite-directed  arcs  of 
capacity  zero.  The  restricted  problem  is  (G',«",c')-  We  call  a  vertex  v,  in  the  restricted  prob¬ 
lem  inner  and  a  vertex  v2  outer ;  we  call  an  arc  (vi  ,v2)  a  split  arc. 

Observe  that  every  circulation  f  in  the  restricted  network  corresponds  to  a  circulation  /  in 
the  target  network  given  by  /(v,  h')=/'(v2,w>j).  By  construction,  the  costs  of  /and/'  are  the 
same.  Conversely,  for  every  circulation  /in  the  target  network,  there  is  a  corresponding  arc  func¬ 
tion/ given  by /'(v2,w1)=/(v,h'), /(v,,v2)=  £  /(x.v),  Function/is  a  circulation  if  and 

(*.*)«  E 
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only  if  /(v  i , v2)  £  m  for  every  split  arc  (v , ,  v2). 

We  assume  (without  loss  of  generality)  that  the  original  problem  has  no  negative  capacities; 
i.e.,  the  identically  zero  arc  function  is  a  circulation.  At  a  high  level,  the  algorithm  consists  of 
initializing  / and  u"  to  be  zero  on  all  arcs  and  repeating  the  following  steps  for  each  bit  of  preci¬ 
sion  in  the  capacities,  proceeding  left-to-right  through  the  bits; 

Step  1.  (Construct  the  new  target  problem  by  introducing  the  next  capacity  bit.)  For  each 
(v.H')  e  £,  replace  u'(v,w)  by  2u'(v,w)  if  the  cunenl  bit  of  «(v,w)  is  zero,  by  2u'(v,h0  +  1  if  the 
current  bit  of  u(v,w)  is  one. 

Step  2.  Construct  the  restricted  problem  (G',u",c'). 

Step  3.  Find  an  optimal  solution  /  to  the  problem  (G'.i/'.c'). 

Step  4.  Construct  the  circulation /in  (G,u',c')  corresponding  to/. 

Step  5.  (Modify  the  target  problem  so  that  the  zero  circulation  is  optimal.)  For  each  (v.w)  e  £, 
replace  u'(v.w)  by  u'(v,w)  -/(v.w). 

The  algorithm  maintains  the  invariant  that  on  entry  to  Step  I,  the  zero  circulation  is  optimal 
for  the  old  target  problem  (as  modified  in  Step  5).  The  following  result  is  similar  to  a  lemma  of 
Gabow  and  Taijan  [8]. 

Lemma  6.1.  The  circulation  /  computed  in  Step  4  is  an  optimal  solution  to  the  target  problem 
(G,u',c'). 

Proof.  By  induction  on  the  number  of  iterations.  On  the  most  recent  entry  to  Step  1,  the  zero  cir¬ 
culation  is  optimal  for  tire  old  target  problem.  This  is  true  for  the  initial  entry  because  of  die  ini¬ 
tialization,  and  true  for  each  subsequent  entry  by  the  induction  hypothesis.  Thus  there  is  some 
price  function  p  for  which  the  reduced  costs  of  all  positive-capacity  arcs  are  non-negative  in  the 
old  target  problem.  By  the  construction  in  Step  1,  all  negative-reduced-cost  arcs  in  the  new  target 
problem  have  capacity  one  or  zero.  Let  /  be  an  optimal  solution  to  the  new  target  problem  such 
that  f  contains  no  zero-cost  cycle  of  positive  flow  (any  such  cycle  can  be  eliminated  by  reducing 
its  flow).  Circulation  f  can  be  decomposed  into  at  most  m  simple  cycles,  each  of  flow  value  one. 
(There  is  at  most  one  such  cycle  for  each  negative  reduced-cost  arc).  The  arc  function  /  on  the 
restricted  network  that  corresponds  to  f  thus  has  /(vj,v2)  S  m  for  each  split  arc  (v,  ,v2),  i.e.  / 
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is  a  circulation.  Since  f  is  optimal  for  the  target  problem,  f  is  optimal  for  the  restricted  prob¬ 
lem.  It  follows  that  any  optimal  solution  to  the  restricted  problem  corresponds  to  an  optimal 
solution  to  the  target  problem,  and  the  lemma  is  true.  □ 

We  shall  describe  an  implementation  of  Step  4  that  uses  the  e-scaling  approach  of  Goldberg 
and  Taijan  [9,11]  (already  discussed  in  Section  3  in  the  context  of  the  transportation  problem) 
with  an  inner  loop  that  is  a  modification  of  the  Ahuja-Oriin  maximum  flow  algorithm  [1,2].  The 
e-scaling  loop  starts  with  a  zero  circulation  and  a  zero  price  function;  the  zero  circulation  is  £- 
optimal  with  respect  to  the  zero  price  function  for  £  =  C.  Then  the  method  iteratively  applies  a 
refinement  subroutine  that  halves  £  and  produces  a  circulation  f  and  a  price  function  p  such  that 

f  is  E-optimal  with  respect  to  p.  When  e  <  the  method  terminates  with  an  optimal  solution 

(Recall  that  the  restricted  network  has  2n  vertices.)  The  formal  definition  of  the  E-scaling  loop  is 
as  follows: 

Step  4.].  Leif  =0.  E  =  C,  and  p-  0. 

Step  42.  For  each  outer  vertex  v2.  let  p(v2)  *  max  {c'(v2,w  j)  I  (v2,w1)e  £'}. 

Step  43.  While  £^cr~.  perform  (c,f.p) «-  refine  (t,f  ,p). 

2  n 

The  special  structure  of  the  restricted  network  allows  the  maintenance  of  the  following 
invariant  in  Step  4.3:  for  any  residual  arc  (v,w)  having  negative  reduced  cost,  v  is  an  inner  vertex. 
We  say  that  a  pseudoflow  / has  the  inner  vertex  property  with  respect  to  a  price  function  p  if  this 
invariant  holds.  Step  4.2  guarantees  that  this  invariant  holds  on  entry  to  Step  4.3. 

The  correctness  of  the  £-scaling  loop  and  the  fact  that  it  terminates  after  OQog(nC))  itera¬ 
tions  of  refine  follow  from  the  results  of  Goldberg  and  Tarjan  [9,1 1]. 

The  following  implementation  of  refine  is  based  on  the  generic  implementation  of  refine 
described  by  Goldberg  and  Taijan  [9,1 1],  specialized  to  use  excess  scaling  as  in  the  Ahuja-Oriin 
maximum  flow  algorithm  [1,2]. 

refine  (tf  ,p). 

Step  R.J.  (Saturate  negative-cost  arcs.)  For  each  split  arc  (vj,v2)  such  that 
c',(v i,v2)<0,  let/(v,,v2)  *  «"(v,.v2). 
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Step  R2.  (Initialize  A,e.)  Let  A  =  2k  where  k  =  f  log  ml .  Replace  e  by  e/2. 

Step  R.3.  (Inner  loop.)  While  A  i  1  repeat  the  following  steps: 

Step  R.3.J.  While  there  is  a  push  or  relabel  operation  that  applies,  perform  such  an 
operation. 

Step  R32.  Replace  A  by  A/2. 

Step  R.4.  (Restore  the  inner  vertex  property.)  For  each  inner  vertex  vj,  replace p(v,)  by 
p(v,)-e/2. 

StepRJ.  Return  (t,f,p). 

The  push  and  relabel  operations  are  defined  as  follows: 
push(v,w). 

Applicability:  ey-(v)  >  A/2,  «/'(w) S  A/2,  u'j-  (v,w)  >  0,  and  Cp(v,w)  <  -e/4. 

Action :  Send  min  {A/2,uy  (v,w)]  units  of  flow  from  v  to  w. 

relabel{\). 

Applicability:  ey  (v)  >  0  and  c'p(v,w)  £  -e/4  for  each  residual  arc  (v.w). 

Action:  Replace  p(v)  by  max  {p(w)-c'(v,w)-e/2). 

Some  remarks  are  in  order  here.  Step  R.1  saturates  all  negative  reduced-cost  arcs,  thereby 
making  f  into  an  optimal  pseudoflow  but  introducing  excesses  and  deficits  at  vertices.  Step  R.3 
moves  the  excess  flow  amounts  to  the  vertices  with  deficits  while  maintaining  e/2-optimality  (for 
the  new  value  of  e).  Step  R.4  restores  the  inner  vertex  property  by  changing  p\  after  this  step,  f 
is  no  longer  c/2-opdmal  but  only  c-optimal. 

Step  4.3  maintains  a  value  A  that  is  an  upper  bound  on  the  largest  excess.  When  no  excess 
exceeds  A/2,  A  is  halved.  All  excesses  are  integers;  by  the  time  A  <  1,  all  deficits  have  been  can¬ 
celed. 

Each  pushing  step  moves  excess  from  a  vertex  with  excess  exceeding  A/2  to  a  vertex  with 
excess  not  exceeding  A/2,  and  through  an  edge  of  cost  between  -e/2  and  -e/4.  Thus  each  push 
is  either  saturating  or  it  moves  at  least  A/2  units  of  flow;  in  the  latter  case  it  reduces  the  cost  of  f 
by  at  least  Ae/8.  We  shall  use  this  cost  reduction  to  bound  the  number  of  nonsaturating  pushes. 
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Wc  have  omitted  a  description  of  how  to  determine  push  and  relabel  operations  that  can  be 
applied.  These  details  can  be  found  in  [2];  the  implementation  uses  the  analogue  to  the  stack- 
pushjrelabel  step  described  there.  The  total  time  spent  in  such  overhead  in  a  single  execution  of 
refine  is  O(nm). 

The  proofs  of  the  following  lemmas  are  easy  modifications  of  proofs  of  analogous  lemmas 
in  [9,10]; 

Lemma  62.  The  push  and  relabel  operations  preserve  e-optimality. 

Lemma  63.  CO  Each  operation  relabel(v)  decreases  p(v)  by  at  least  e/4.  00  During  an  execu¬ 
tion  of  refine ,  the  maximum  amount  by  which  p(y)  can  decrease  is  6tn.  for  every  vertex  v.  (iii) 
The  total  number  of  relabel  operations  during  an  execution  of  refine  is  0(n2),  taking  0(nm ) 
time.  0V)  The  number  of  saturating  pushes  during  an  execution  of  refine  is  O(nm). 

Lemma  6.3  imphes  that  the  running  time  of  refine  is  0(nm )  plus  (9(1)  per  nonsaturating 
push.  We  shall  establish  a  bound  of  0(n2  log  m  +  nm)  on  the  number  of  nonsaturaling  pushes, 
thereby  obtaining  an  0((n2  log m  +  nm)  log  V  log (nC))  bound  on  the  triple  scaling  algorithm. 

To  bound  the  number  of  nonsaturaling  pushes,  we  define  the  cost  of  a  pseudofiow  f  with 
respect  to  a  price  function  p  to  be 


costp(f ')  —  X  cp(v,w)f(v,w)  =  £  c'(v,w)f(v,w)~  £  ef(v)p(v). 

/'(»,»)  >  0  /'(»,<•>)  >o  *  «  v 


Observe  that  it  f  is  a  circulation  its  cost  does  not  depend  on  p ,  i.e.  cosXplf’)  =  cost( /')  where 

costif )«  ]£  c/(v,w)  f* (v.h')  (as  defined  in  Section  2 ). 

/'(».»)  >o 

As  noted  above,  a  nonsaturating  push  decreases  costp(f)  by  at  least  Ae/8.  A  relabeling  of  a 
vertex  v  that  decreases  ply)  by  an  amount  x  increases  cosiplf)  by  zey(v). 

We  want  to  arudy  how  much  costplf)  can  vary  during  an  execution  of  refine.  We  do  this 
by  relating  cost plf)  to  the  cost  of  an  optimal  circulation.  This  requires  the  following  lemma, 
which  states  a  general  result  about  circulations  and  pseudoflows. 


Lemma  6.4.  Let  /  be  a  circulation,  f  a  pseudoflow,  and  p  a  price  function  cm  a  network  with 
capacity  and  cost  functions  u  and  c,  respectively.  Then 
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COSlpif  )  -  COSt(f)  £  £  -cr(v,w)  iy(v,w). 

(».W).  Cj(y.w)  <  0 


Proof  Let  f  be  the  pseudoflow  obtained  from  f  by  saturating  aU  negative  reduced-cost  ares. 
Then 


costpif)  s*  costpfj*')  •+  £  -Cp(v,w)Uf’(v,w). 

(»,»)  «,(*.  ')<  o 

But  /can  be  obtained  from  f  by  increasing  flow  along  a  collection  of  paths  and  cycles  of  arcs  in 
C/„,  each  of  which  has  nonnegative  reduced  cost  Thus  cost(f)  £  costp(f"),  and  the  lemma  fel¬ 
lows.  □ 

Consider  a  time  during  the  execution  of  refine.  For  the  current  pseudoflow  f  and  price 
function  p ,  we  define  a  potential  <S>  by  <X>  =  (costp(/')  -  cost  (f ))  /(Ae),  where  f  is  any  optimal 
circulation.  The  following  lemma  bounding  O  is  the  heart  of  the  analysis  of  nonsaturating 
pushes. 

Lemma  65.  — 1 t&n2  S  <J>  S  6nm/A. 

Proof.  We  shaJ  prove  that  -48Ae/i2  £  cost^Cf)  -  cost(/*)  S  6enm.  To  obtain  the  upper  bound, 
we  note  first  that  f  is  e/2-optimal  with  respea  to  p.  By  decreasing  prices  on  inner  vertices  by 
e/2,  we  obtain  a  price  function  p'  with  respea  to  which  f  is  e-optimal  and  such  that  f  has  the 
inner  vertex  property  with  respect  to  p' .  Since  no  vertex  has  excess  exceeding  A  S  2m,  cosy  C/*) 
icostp(/')  +  wm.  By  Lemma  6.4, 


cost p’(J')-cost(J*)£  J)  -c'p’  (y,w)Uf'(v,w). 

(»,w>  <  o 


The  sum  of  aO  positive  excesses  is  at  most  2/iA  S  4 wr,  hence  ao  is  the  negative  of  the  sum 
of  all  negative  excesses.  For  any  inner  vertex  v,  the  total  residua]  capacity  of  all  arcs  leaving  is 
m-e)(v).  Hence  the  total  residual  capacity  of  all  arcs  leaving  inne*-  vertices  is  at  most 
nm  +  4am  =  Snm.  By  the  inner  vertex  property,  every  negative  reduced-cost  arc  leaves  an  innei 
vertex.  Hence 
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£  -Cp’(y,w)uf  (v.w)  £  5enm. 

(»,»•):  c''  (r.w)  <  0 


Combining  inequalities  gives  costp  if)-  cost(f)H.  6 tnm. 

To  obtain  the  lower  bound,  we  note  first  that  we  can  assume  that  there  is  a  price  function  p’ 
with  respect  to  which  f  is  optimal  arid  such  that  lp(v)  -p*(v)  I  <  12 tn  for  every  vertex  v.  This  is 
since  the  repeated  executions  of  refine  in  Step  4.3  of  the  e-scaling  method  will  produce  an 
optima]  circulation  f  and  a  final  price  function  p*  such  that  f  is  optimal  with  respect  to  p*;  the 

m 

total  price  change  of  any  vertex  over  all  the  iterations  of  refine  is  at  most  £  6nt/2‘  *  \2cn  We 

iM) 

have  co$t{f*)-costp-{f)£  0,  since  f  can  be  converted  into  f  by  increasing  flow  on  residual 
arcs  of  Gy.  all  of  which  have  nonnegative  reduced  cost.  But 

costp-(f)-costp(f)£  J  lp*(v)-p(v)lley(v)l£48AEn2 

*tv 


since  the  sum  of  the  positive  excesses  is  at  most  2nA,  as  is  the  negative  of  the  sum  of  the  negative 
excesses  Combining  inequalities  gives  cost{f )  -  costp(f )  £  48Ae/i2.  □ 

Lemma  6.7.  The  number  of  nonsaturating  pushes  during  an  execution  of  refine  is 
G(n2  logm  +  nm). 

Proof.  Each  nonsaturating  push  decreases  C>  by  at  least  1/8.  Any  saturating  push  also  decreases 
<I>.  A  relabeling  that  decreases  the  price  of  a  vertex  by  x  increases  $  by  at  most  x. 

There  are  at  most  log  m  +  1  iterations  of  Step  R.3  in  refine.  Consider  the  i,k  iteration.  Sup¬ 
pose  that  during  this  iteration  the  total  decrease  in  vertex  prices  is  X;.  If  p,  is  the  number  of  non¬ 
saturating  pushes  during  this  iteration,  then  the  iteration  causes  a  net  decrease  in  C>  of  at  least 
p,78-x Changing  A  between  two  iterations  can  increase  C>  by  at  most  a2.  Summing  over  all 
phases  ard  applying  Lemma  6.5,  we  have 


X(Pi/8-x,)S6nm  +48n2  +  n20ogm+  1). 


Since  £x,  *  0(n 2),  it  follows  that  £p,  *  0(n2  log  m  +  nm).  □ 

i  i 


Theorem  6.8.  The  triple  scaling  algorithm  for  finding  a  minimum-cost  circulation  runs  in 
0((n2  log  m  +  run)  logt/  log(nC))  time. 

Using  dynamic  trees  in  the  inner  loop  of  the  triple  scaling  algorithm,  analogously  to  the 

way  they  are  used  in  [2],  reduces  the  time  bound  to  0(nm  log  (1  +  —  logm)  log U  log  (nC)). 

rn 

Further  miqpr  improvements  might  be  possible  using  additional  ideas  in  [2J.  We  shall  not  pursue 
this  possibility  further,  however,  since  in  any  case  the  approach  of  Sections  3-5  produces  better 
bounds. 

We  conclude  this  section  by  noting  that  there  is  an  alternative  way  to  solve  the  sequence  of 
restricted  problems  generated  in  the  outer  capacity-scaling  loop.  Namely,  we  can  use  a  simple 
version  of  the  network  simplex  rule,  specifically  Damzig’s  minimum  reduced-cost  pivot  rule  with 
lexicography  to  avoid  cycling.  This  rule  was  studied  by  Oriin  [15],  who  obtained  a  bound  of 
0(nmU  log  (mUC))  on  the  number  of  pivot  steps.  In  a  restricted  problem,  U  =  m,  and  the 
capacity-scaling  method  with  Step  3  implemented  using  the  network  simplex  method  runs  in 
0(nm  3  log  U  log  (mC))  time,  if  the  time  to  do  one  pivot  step  is  0(m).  In  this  algorithm,  it  is  not 
necessary  to  transform  the  graph  by  splitting  vertices;  it  suffices  to  impose  a  capacity  bound  of  m 
on  every  arc  in  G.  This  time  bound  is  not  noteworthy;  more  interesting  is  the  mere  fact  that  com¬ 
bining  one  scaling  loop  with  a  standard  version  of  the  network  simplex  algorithm  gives  a 
polynomial-time  algorithm.  Without  the  scaling  loop,  the  same  version  of  the  network  simplex 
algorithm  car  take  exponentially  many  pivot  steps  [24]. 

7.  Remarks 

We  have  shown  that  the  minimum-cost  circulation  problem  can  be  solved  in 
0(nm  log  log  U  log  (nC))  time,  and  even  in  0(nm  log  (2  +  —  log  l/)  log  (nC))  time  if 

fit 

log  G  *  n  (—  log  n).  We  have  derived  analogous  bounds  for  the  transportation  problem.  Our 
n 

algorithms  use  scaling  of  both  costs  and  capacities,  combined  with  an  augmenting  path  method 
and  an  implementation  based  on  dynamic  trees.  If  dynamic  trees  are  not  used,  the  time  bound  for 
the  minimum-cost  circulation  problem  is  0(nm  log  C/ (I  +  log  (nC)/ log  log  f/)).  or 

O  (n2  Jog  V  (1  ♦  log  (nC)/  log  log  17))  if  log  U  =  O  (—  logn).  Under  the  similarity  assumption 

It 

(7],  namely  \o%U  *0(k>gn)  and  log  C  ■  0(logn),  our  time  bound  with  dynamic  trees  is 

„2 

0(nm  logn  loglogn),  which  beats  the  best  previous  bound  of  0(nm  log n  log( — ))  [11]  except 

m 

n2 

for  very  dense  graphs  (m  *  Q  (- - )). 

logn 
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We  expect  that  some  version  of  our  algorithm,  probably  without  dynamic  trees,  will  in  prac¬ 
tice  be  competitive  with  or  superior  to  previously  existing  algorithms.  We  have  not  yet  done  the 
required  experiments  to  confirm  or  refute  this  hypothesis.  Our  experiments  with  similar  algo¬ 
rithms  (see  e.g.  [9])  suggest  that  periodic  scans  to  tighten  prices  may  increase  the  practical, 
though  not  the  theoretical,  speed  of  the  algorithm. 

We  have  discussed  a  capacity-bounding  technique,  which  allows  us  to  use  a  modification  of 
the  Ahuja-Orlin  maximum  flow  algorithm  in  the  inner  loop  of  the  Goldberg-Tarjan  minimum- 
cost  flow  method.  The  analysis  of  this  technique  uses  the  cost  of  the  current  pseudoflow  as  a 
measure  of  its  quality.  This  makes  the  analysis  very  intuitive. 

A  tantalizing  open  question  is  whether  there  is  an  O  ( nm  log  logn  log  (nC))-time  algorithm 
for  the  minimum-cost  circulation  problem.  We  believe  that  the  answer  is  yes  and  that  a 
modification  of  our  methods  will  lead  to  such  a  bound.  Such  a  result  would  probably  give  a  time 
bound  of  O  (nm  log  logn)  for  the  maximum  flow  problem,  which  would  also  be  an  improvement 
over  known  results.  (See  [2,10].) 
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